TRUNCATED MOMENT PROBLEMS IN THE CLASS OF 
GENERALIZED NEVANLINNA FUNCTIONS 
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Abstract. Truncated moment problems in the class of generalized Nevan- 
linna functions are investigated. General solvability criteria will be established, 
covering both the even and odd problems, including complete parametrizations 
of solutions. The main new results concern the case where the corresponding 
Hankel matrix of moments is degenerate. One of the new effects which reveals 
in the indefinite case is that the degenerated moment problem may have in- 
finitely many solutions. However, with a careful application of an indefinite 
analogue of a step-by-step Schur algorithm a complete description of the set 
of solutions will be obtained. 



1. Introduction 

The main purpose of this paper is to study general truncated (real) moment 
problems and some associated interpolation problems involving a finite sequence of 
real numbers So, Si, . . . , si. In order to describe some of the contents and results in 
the paper it is natural to start by recalling a couple of notions and results appearing 
in classical truncated moment problems. 

The truncated Hamburger moment problem for real numbers So, s ij • ■ • j &l € 
Z_|_) consists of finding a positive measure fi on / = R for which 



(1.1) Ji?dn(t)=8i, j = 0, 

This problem will be called odd or even , if the number I is odd or even, respectively. 
In the case where / = R+ the problem (jl.lj) is called the truncated Stieltjes moment 
problem. Due to the Hamburger- Nevanlinna theorem [TJ Theorem 3.2.1] in the even 
case (£ = 2n) the conditions in (jl.lj) can be rewritten in terms of the associated 
function <p(\) defined by 

(1.2) p(A)= / ^% A e C \ R, 



ht-X' 

as the following interpolation problem at A = oo: 

(1-3) <p(X) = ~ WT + °{wi)' A ^°°' 

The notation A^oo means that A — > oo nontangentially, i.e. S < arg A < tt — S for 
some 5 > 0. Recall, that (p belongs to the class Nq of Nevanlinna functions, i.e., 
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ip(X) is holomorphic on C \ R, satisfies the symmetry condition ip(\) = </?(A), and 
has a nonnegative imaginary part for all A G C+. The moment problem can 
now be reformulated as follows: find a Nevanlinna function <f(X) such that (|1.3[) 
holds. It follows easily from (jl.ip that the following inequality 

(1-4) S n := (s i+j )l j=0 > 0, 

is necessary for the problem to be solvable. In the case where the matrix 

S n is invertible this condition is also sufficient for to be solvable, and all its 
solutions are described by the formula (see [28]): 

n ^ f Mt) _ Q»(A)r(A) + Q n+1 (A) 

1 ' j ^ [ ' J R t-\ P n (X)r(X) + P n+1 (X) ' 

where P„ are polynomials of the first kind orthonormal with respect to (3, 

'P„(f)-P n (A) > 



(1-6) Q»(A) = 6 , / _ v 

are polynomials of the second kind and r(A) is an arbitrary Nevanlinna function 
from the class No which satisfies the Nevanlinna condition 

(E) t(A) = o(A) as A^Foo; 
in (|1.6[) 6 stands for the nonnegative functional defined on the set C[t] of polyno- 
mials via 

G@) = 8j, j = 0,l,...,2n. 
Notice that this classical result depends essentially also on the assumption that 
the moment problem is even (i.e., I = In). Indeed, the odd Hamburger moment 
problem is not equivalent to the interpolation problem (jl.3l) . A convenient frame- 
work to formulate the problem in the odd case is provided by the classes No,-^ 
appearing in [21]: they consist of functions / G No of the form (|1.2j) . such that the 
measure u in (|1.2j) satisfies the condition 

(1 + \t\ e )da(t) <oo (£= -1, 0,1,2,... ). 

Then the Hamburger-Nevanlinna theorem can be restated as follows: a is a solution 
of the moment problem (|l.ll) if and only if the associated function ip belongs to 
the class No _£ and has the asymptotic expansion (| 1 . 3[) . It is a consequence of 
the results in the present paper (see Corollary I5.2[) that the set of solutions of the 
nondegenerate odd moment problem (ll.l[) can also be given in the form (jl.5|) , where 
r now ranges over the class N 0i i and satisfies 

(O) r(A) = o(l) as A^oo. 
It should be mentioned that in the special case of the nondegenerate odd Sticltjcs 
moment problem the set of solutions was parametrized in [23 with the parameter 
r ranging over the class of Stieltjes functions (i.e. Nevanlinna functions of the 
form (II. 2[) with / = K+). Such functions automatically belong to the class N .i 
and they satisfy also the condition (O). However, it seems to the authors that 
for general measures, whose support is not contained in some semiaxis in R, the 
above mentioned description of the solution set for nondegenerate odd Hamburger 
moment problem has not appeared in the literature earlier. 

In the case where the matrix S n is degenerate the condition (|1.4j) is not anymore 
sufficient for the problem (|1.1[) to be solvable; see [26, 22, 8 . Recall that the Hankel 
rank, denoted by rank (s, 2n), of the sequence (s, 2n) = {s^}?™ is defined as follows: 
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rank(s, 2n) = n+1 if detS n ^ 0, otherwise, rank(s,2n) is the smallest integer r, 
< r < n, such that 



/ s r 

(1.7) : 



€ span 



V 



/ S r - 



\S r -l 




By a Frobenius theorem (see [19j Lemma X.10.1]) Hankel rank of (s, 2n) is the 
smallest integer r, 1 < r < n + 1, such that 

det 5 r _i 0, and det Sj = for j > r. 

In particular, rank(s, 2n) = if sq = ■ ■ ■ = s n = 0, otherwise rank(s,2n) is the 
smallest integer r, 1 < r < n + 1, such that 

det S r -i 7^ 0, and det Sj — for j > 7-. 

A sequence (s,2n) = {sj}|™ with the Hankel rank r — rank (s, 2n) is called 
recursively generated, if there exist numbers ao, . . . , a r -i, such that 

(1.8) Sj = a Sj^ r H + a r _iSj-_i (r < j < 2n). 

Now solvability criteria for the degenerate truncated moment (|1.1[) can be for- 
mulated as follows. 

Theorem 1.1 ( 8 ]). Let the matrix S n — (sj+j)™ J=0 &e nonnegative and degenerate, 
and let r = ranks. Then the following statements are eguivalent: 

(i) the moment problem (|1.1[) is solvable; 
iii) rankS'n = r; 

(hi) S n admits a nonnegative Hankel extension S n +i; 

(iv) the sequence (s,2n) — {sj}|™ is recursively generated. 

If any of the assumptions (i)-(iv) are satisfied, then the problem (jl.lll has a unique 
solution (p(\) = — p r p^ ■ 

It is interesting to note that Theorem 11.11 contains as a corollary the following 
rigidity result due to D. Burns and S. Krantz [7]: if y is a rational Nevanlinna 
function of degree r with the asymptotic expansion (|1.3I) and ip is a Nevanlinna 
function such that tp(X) = <p(\) + o( A ar+i ) as A^oo, then -0(A) = ^(A). 

The main subject of the present paper is the study of degenerate odd and even 
moment problems involving hnite sequences Sq, s±, . . . , si of real numbers by means 
of functions belonging to the class of generalized Nevanlinna functions, which con- 
tains the class of Nevanlinna functions appearing in (jl.2p as a subclass. 

Definition 1.2. ([24]) Let k E N. A function (p meromorphic on C+ is said to 
be from the class N K , k, € N, of generalized Nevanlinna functions with k negative 
squares, if the kernel 

MA) = 

A — uj 

has k negative squares on C+, i.e. for every choice of m € N, Ai, . . . , A m € C+ the 
matrix 

(NA fc (A,))r fe =i 

has at most k and for some choice of m, Xj exactly k negative eigenvalues. 
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In [T2J [14] (see Definitions 12. 31 l2~7l below) subclasses N Kj _^ of the class N K were 
introduced as indefinite analogues of the subclasses No,-^ appearing in [21] . 

In this paper we consider in a parallel way the following two problems: 
Indefinite truncated moment problem MP K (s,£): Given are k € N, £ £ Z+, 
and So, . . . , Si 6t. Find a function <p G N K ._£ with the asymptotic expansion (jl.3l) . 
Denote by M. K (s,£) the set of solutions of this problem. 

Multiple indefinite interpolation problem IP K (s, £): Given are k £ N, £ £ Z + , 

and so,..., st € K. Find a function ip € N K with the asymptotic expansion (| 1 . 3|) . 
The set of functions with these properties is denoted by T K (s,£). 

As was mentioned above, the problem MPq(s,1) is equivalent to the trun- 
cated Hamburger moment problem Furthermore, in the even case (£ = 2n) 
M K (s,£) — T K {s,i), while in the odd case (£ = 2n+l) we have M K (s,£) C l R (s,£), 
but the reverse inclusion fails to hold in general. 

The problems MP K (s,£), IP K (s,£) will be called nondegenerate, if 

(1.9) detS„ ^ 0, for n = [1/2]; 

otherwise they are called degenerate. The indices j for which detSj_i ^ are 
called normal indices of the Hankel matrix S n . Let 

n\ < ri2 < • • • < 7t.jv < n + 1 

be the sequence of all normal indices of the matrix S n . In the case of arbitrary 
Hankel matrix S n we show that the largest normal index njy of S n coincides with 
the Hankel rank of the sequence s. 

A necessary condition for the problems MP K (s,£), IP K (s,£) to be solvable is 
that 

(1.10) K>V-(S n ), 

where v-{S n ) is the total multiplicity of all negative eigenvalues of S n . The method 
we use for the solution of the moment problem MP K (s,2n) and IP K (s,£) for 
K > U-(S n ) is based on the Schur-Chebyshev recursion algorithm, studied in the 
nondegenerate situation by M. Derevjagin [9] (see also [4]). With this method every 
solution ip of the moment problem MP K (s, 2n) can be obtained via 

(1.11) tp(X) 



pi(\) + ei<pi(\) ' 

where pi (A) = P ni (X), £\ = sgns ni _i, and ip± is a solution of an "induced" moment 
problem MF K _ Sl (sW,2(n- n\)) with k\ — V-(S ni -i). 

In the case of a nondegenerate moment problem the condition (| 1 . 10[) is also 
sufficient for the problem MP K (s,£) to be solvable and subsequent applications of 
the formula (|1.11[) shows that M K (s,l) in the even case (£ = 2n) is parametrized 
via the linear fractional transformation 

(A)r(A)+Q njv (A) 

(L12) ^^iOWA)TiUA)' 

with the parameter r ranging over the class ^ K -v^(s n ) an d satisfying the Nevan- 
linna condition (E) ; see [T5] [TT] . In this formula Pj and Qj are polynomials of the 
first and the second types introduced in [TT]. In the odd case a similar description 
of the sets .A/f K (s,£) and I K (s,£) is given in Theorem I5.ll with the parameter r 
ranging over the class ^ K - V _(s n ),i and N K _ 1/ _(5 ji ), respectively, and satisfying the 
condition (O). It should be mentioned, that this result for X K (s, £) can be derived 
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also from the recent paper [5] on boundary interpolation in generalized Nevanlinna 
classes. 

Now let us briefly describe some of the main results obtained in the present paper 
for the degenerate indefinite truncated moment problem in the even case. As in the 
definite case, for a degenerate problem the condition (II . 10[) is not sufficient for the 
problem MP K (s,2n) to be solvable. The following theorem gives some solvability 
criteria in the special case, where k = v-(S n ); in fact, this result offers a natural 
generalization for the results due to Curto and Fialkow [5] , which were formulated 
in Theorem II. f I above. 

Theorem 1.3. Let ri\ < ri2 < ■ ■ ■ < be the sequence of all normal indices 
of a degenerate matrix S n — (si+j)™ J=0 and let k — V-(S n ). Then the following 
statements are equivalent: 

(i) the moment problem MP K (s,2n) is solvable; 

(ii) rankSVj = njv; 

(iii) S n admits a Hankel extension SVi+i such that 

V-(S n+1 ) = V-(S n ); 

(iv) the sequence s = {sj}^™ is recursively generated. 

If any of the assumptions (i)-(iv) is satisfied, then the problem M P K (s, 2n) has a 
unique solution <p{X) = — p""^) ■ 



As a consequence of Theorem 11.31 a rigidity result for generalized Nevanlinna 
functions from [5] can be derived (see also [5]). 

A new effect which appears in the indefinite case is that the degenerate moment 
problem MP K (s,2n) has infinitely many solutions for k large enough. As will be 
shown below, the problem MP K (s, 2n) with k > v^(S n ) is solvable if and only if 

(1.13) k >v-{S n ) + v {S n ). 

If, in addition, ranks'.,! = rijv + 1 and v satisfies some appropriate further condi- 
tions (see (|2.17l) below), then the solution set A4 K (s, 2n) can be described by the 
formula (|1.12[) . where 

< L "> «*> = % 

and t is a function from the class N K _„, which satisfies (E); see Theorem 15.81 

On the other hand, if (|1.13[) holds and rankSVi > + 1, then the solution set 
M K (s, 2n) is described by (|1 . 12|) . where 

(L15) T(A) = A^(g(A)+ e r(A)) ' 

p is a polynomial of degree n — vq + 1 (as given in (j3. 15[) below) , and r is a function 
from the class N^.,,, which satisfies (E); see Theorem 15.91 

In the odd case the degenerate indefinite truncated moment problem can be 
treated analogously. A condition, similar to (|1.13[) appeared in [55] as a solvability 
condition for a degenerate indefinite Nevanlinna-Pick interpolation problem. 

In Section 2 the basic tools needed in this paper are given. Solutions to the 
so-called basic moment and interpolation problems will be described in Section 3. 
Section 4 describes a general Schur-Chebyshev recursion algorithm, which makes 
use of the normal indices of the associated Hankel matrix S n = (sj+i)?j=o defined 
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in (|1.4|) . In Section 5 solvability criteria and complete descriptions for the set 
of solutions of the problems MP K (s,£) and IP K (s,£) in the general setting are 
established. Finally, Appendix contains some results on block matrices, which 
are needed in this paper; however, they may be also of independent interest: for 
instance, Lemma lA.21 gives an extension of a well-know result on nonnegative block 
matrices. 

2. Preliminaries 

2.1. Canonical factorizations of generalized Nevanlinna functions. The 

definition of the class N K of generalized Nevanlinna functions is given in the Intro- 
duction. Clearly, if ip G N K and <p ^ 0, then also —l/ip£ N K . 

Recall (see |25j). that the point a G R is called a generalized pole of nonpositive 
type (GPNT) of the function ip G N K with multiplicity K a (p) if 

(2.1) - oo < lim(z - a) 2Ka+1 ip(z) < 0, < lim (z - af Ka ~ l ip(z) < oo. 

z — >a z— >a 

Similarly, the point oo is called a generalized pole of nonpositive type of tp with 
multiplicity Kqo(p) if 

(2.2) < lim < oo, -oo < lim ?( z ) < 

A point /? 6 RU {oo} is called a generalized zero of nonpositive type (GZNT) of 
the function cp G N K if (3 is a generalized pole of nonpositive type of the function 
—1/ip. The multiplicity Kf}{<p) of the generalized zero of nonpositive type /3 of tp 
can be characterized by the inequalities: 

rt - n . v Viz) ^ / v ^ n 

(2.3) < lim — : < oo, -oo < lim — r- — T < 0. 

V ; z-=f/3 (z - l3) 27T f + 1 - z^f3 (z - /S) 2 ^- 1 - 

Similarly, the point oo is a generalized zero of nonpositive type of tp with multiplicity 

7Too(</>) if 

(2.4) - oo < lim z 2 *°° +1 tp{z) < 0, < lim z 2lx ^- x p(z) < oo. 

z— foo z— >oo 

Remark 2.1. If ipi G N Kl and p>2 G N K2 then ipi + ip2 belongs to N K , where 
K < Ki + K2- It was shown by M.G. Krem and H. Langer in [25] for <p G N K that the 
total multiplicity of poles (zeros) in C+ and generalized poles (zeros) of nonpositive 
type in RU {oo} is equal to n. As a corollary of this result one obtains that if 
p>\ G N Kl and ip2 G N K2 have no common poles in C+ and common generalized 
poles of nonpositive type in lU {oo} then ipi + p>2 G N Kl + K2 . 

The generalized poles and zeros of nonpositive type of a generalized Nevanlinna 
function give rise to the following factorization result ([IT], see a ls° [12j). 

Theorem 2.2. Let tp G N K and let a%, . . . , ai (7?i, ■ ■ ■ , /3 m ) be all the generalized 
poles (zeros) of nonpositive type of 'ip inM. and the poles (zeros) of <p in C + with mul- 
tiplicities Ki, . . . , At; (Vi, . . . , TT m ). Then the function p admits a (unique) canonical 
factorization of the form 

(2.5) p{z) = r(z)r*(z)p {z), 

where tpo G No, r#(z) = r(z), and r = p/q with relatively prime polynomials 

rn I 

POO = T[(z - PjY*, ?(*) = ]J(z - aj)"', 
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of degree n — TToo(ip) and K — Koo(</9) 7 respectively. 

It follows from (|2.5|) that the function ip admits the (factorized) integral repre- 



sentation 



(2.6) ,( Z )=r( Z )rnz)(a + bz + j^ z - Y ± w ) 



P 

d P(t) ) > r = -) 
9 



where a e 1, i > 0, and p(£) is a nondecreasing function satisfying the integrability 
condition 

dp(t) 
t 2 + l 



( 2 - 7 ) / <°°- 



2.2. The subclasses N re ,_^ of generalized Nevanlinna functions. 

Definition 2.3. (see [T2]) A function £ N K is said to belong to the subclass 
X, .,. h 

(2.8) lim *M = and f ' Im ^ ' d y < oo, 

z^oo 2 J ?) y 

with 77 > large enough. Similarly, a function 6 N K is said to belong to the 
subclass N Ki0 , if 

(2.9) hm ^ - — and lim sup |zlm<^(z)| < 00. 

z— >oo Z ^^Tqq 

Remark 2.4. Every function p 6 N Kj i has a nontangential limit lim,^^ <p(\) at 
infinity. As was shown in |12j the following implication holds: 

^eNu, lim <p(\) ^ =► G N M . 

A— >oo 

In the following theorem the subclasses N Kj i and N K .o are characterized in terms 
of the integral representation ()2.6[) . 

Theorem 2.5. (|12j) .For ip £ N K and £ = 0, 1 the following statements are equiv- 
alent: 

(i) p belongs to N K j; 

(ii) 93 /ias the integral representation (|2.6[) wii/i degg — degp = ^^((p) > 7 or 
with degp = degg (it a3 { l P) = ^)> b — 0, and 

(2.10) f (1 + \t\)- e dp(t) <oo. 

Jr 

Remark 2.6. If tp € N Ki o then the statement (ii) in Theorem l2.5l can be strengthened 
in the sense that for every function tp € N Ki o there are real numbers 7 and sq, such 
that 

(2.11) (p( z )= 7 - ?± + (Vj , z ^oo. 

Definition 2.7. (Q3]) A function 95 € N K is said to belong to the subclass N K ,_2n, 
n 6 N, if there are real numbers 7 and So, . . . , S2 n — 1 such that the function 



(2.12) ^ 2 "J(z) = z 2 "UM-7 + E 



2n 



is 0(1/ z) as z^oo. Moreover, 93 € N K is said to belong to the subclass N K) _2n+i 
if the function ip^ in (|2.12p belongs to N«/,i for some «/ € N, «' < k. 
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As was shown in [14], the following inclusions are satisfied 

(2.13) • • • C N K ,_ 2 „-1 C N«,_2„ C N« _2n+l C • • • C N K , C N M . 

The subclasses N K) _.g, I £ N, can be characterized by means of the integral 
representation of p in (12 .6[) . 

Theorem 2.8. /|14j ) i^or ip £ N K £/ie following statements are equivalent: 

(i) p £ N K ,_ £; £ e N; 

(ii) (/? /ias an integral representation (|2.6[) tuif/i Tr^ty?) = degg — degp > (and 
b = if TToc{(p) — 0), such that 

(2.14) / (i + |;|) £ - 2 ^a>(;) <oo. 



Remark 2.9. By Definition 12.71 every function 9? £ N Ki _ £ with odd ^ admits the 
asymptotic expansion 

(2-15) ^(A)= 7 -| _g_ +0 A^oo. 



If <p £ N K ._£ and £ is even due to Theorem 12.81 there exists a real number s 2 n, 
such that (|2.15[) holds. Conversely, if <p £ N K and satisfies (|2.15[) for even £, then 
p £ N Kj _^. This proves that in the even case (I — 2n) SA K (s,£) = X K (s,£), while 
in the odd case (£ = 2n + 1) the M K (s,£) C l K (s,£). 

The following Lemma is immediate from Definition 12. 71 Theorem 12.81 and Re- 
mark mu 



Lemma 2.10. Let ip £ N K; _^ let Vq < min{ir oa (ip),£/2} and let 
(2.16) 0(X) = X 2v °p{X). 

Then ip £ N K _ v _^_ 2vo ) , where 

, 2 17) v= { V > °! 

Conversely, if £ N- and Koo(p) = 0, ip and v are given by (j2 . 16[) and (|2.17p . 
i/ien <y3 e N~ +j/ __ ( f +2j/o) and ^0 < ttoc^)- 

Proof. 1) Let y> e N K) _^. Due to Theorem 12.81 ip admits the canonical factoriza- 
tion (|2.5j) . where the measure p in the integral representation (|2.6|) satisfies the 
condition (|2.14p . It follows from (|2.5I) and (|2.17l) that admits the canonical fac- 
torization 

p(A)p#(A) 

^ = PF(A)" 0() ' 
with the same function pq £ No and 

TToo(p) = degg - degp = ^(v?) - 
Hence, the condition (|2 . 14|) takes the form 

(2.18) / {l + \t\f- 2 ^\dp{t) <oo. 



where £ = £ — 2vq. Due to Theorem 12.81 the latter condition is equivalent to the 
inclusion 

p £ N Ki _^_ 2l/0 ) 
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with k < K. 

2) Since vq < Hoo(p) then neither p nor p has a GPNT at oo. Assume that 
Ko{ip) > vq. Then both tp and tp have GPNTs at and Ko(p) = kq(p) + vq. 
Counting the total pole multiplicities of p one obtains tp £ N K _„ by Remark 1 2. 11 

Assume that < Ko(tp) < vq. Then tp has no GPNT at and tp has a GPNT at 
of multiplicity kq (tp). Therefore, <p € N K _ K[) ( ¥ ,). 

And finally, if kq(p) = and ttq(p) > 0, then neither tp nor tp have a GPNT at 
and thus tp € N K . All the above statements are easily reversed. □ 

2.3. Toeplitz matrices. A sequence (c, n) :— (co,...,c„) of (real or complex) 
numbers determines an upper triangular Toeplitz matrix T(cq, . . . ,c n ) of order 
(n + 1) x (n + 1) with entries t%j = Cj^i for i < j and tjj = for i > j. 



V 



(2.19) T(c m ,..., Cj ) 
Clearly, T(c m , . . . ,Cj)* = Jj- m+ xT(c m , 

(2.20) Jj-m+l = 



'0 V 



< m < j < n. 
-i, where 



1 Oy 

In particular, the coefficients of the asymptotic expansions 



(2.21) 



c(A) = c + 
d{\) = d + 











A n 










" + A7~ 



determine the Toeplitz matrices 
/ c 

(2.22) T(c ,... ) c n ) = 

V 



A— >oo; 
A^>oo. 

/ do 



T(d , ...,d n ) 



Co 



V 



d 



Lemma 2.11. Le£ i/ie functions c and d (meromorphic on C\M.) have the asymp- 
totic expansions (|2.2ip and /ei a(A) = c(A)d(A) have the asymptotic expansion 



a(A) = a o + — H h A" +0 



A' 



A^oo. 



Then the first n + 1 coefficients of the asymptotic expansion of a(A) can be found 
by 

(2.23) T(a , . . . , a„) = T(c , . . . , c„)T(d , . . . , d„). 

Lemma 2.12. The formula 



(2.24) 



P(A) 



det 5„ 



det 



V i 



A 



S2m+1 

A m+i y 
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where Sj = 0,j < m, s m ^ (m > 0), and S m is as in (|1.4[) . defines a monic 
polynomial p(X) — 5Zy^g PjA^ of degree m+1 whose coefficients pj , 1 < j < m + 1, 
satisfy the matrix equality 

(2.25) T(p m+ i, . . . ,px)T(s m , . . . , S2m) = s m I m +i (Pm+l = 1) 

for an arbitrary real number S2 m +i- 

Proof. Evaluating the determinant in (|2.24[) with respect to the last row shows im- 
mediately that p(X) is a monic polynomial of degree m+1. To see that the coeffi- 
cients pj of p(A) in (|2.24l) satisfy (|2.25[) substitute A-? by Sj + k-i for j = 0,1,..., m+1 
(s_i = 0) in the formula (|2.24p . Then for k = the evaluation of the determinant 
in (|2.24p yields the equality 

V m+1 Ti « , — a 

and for k = 1, . . . , m one obtains 

j=i ^Sj+fe-i = 0. 

This means that the polynomial p defined by (|2.24|) automatically satisfies (|2.25|) 
for arbitrary S2m+i S K. Note that S2m+i only appears in the constant coefficient 
Po of p(X), which can be seen e.g. by evaluating the determinant in (12.24[) with 
respect to the last column. □ 

2.4. Asymptotic expansions of certain fractional transforms. In the next 
lemma the polynomial p(X) defined in Lemma [2.121 (see ()2.24|) ) appears when in- 
verting an associated asymptotic expansion. 

Lemma 2.13. Let (s,£) = (sj)j =0 be a sequence of real numbers such that Sj = 0, 
j < m, and s m ^ 0, £ > 2m, and let the monic polynomial p(X) — ^2™^ PjX-' be 
defined by (|2.24[) . Then a function ip (meromorphic on C\RJ admits the asymptotic 
expansion 

(2-26) ^--^-•••-A^ +0 (a^)' 

if and only if the function —s m /(p(X) admits the asymptotic expansion 

(2.27) - s m /<p(X) = P (X) + er(A) A^>oo, 

where e = sgn s m and t(X) satisfies one of the following conditions: 

(i) if £ = 2m then t(A) = o(X), X^oo, and in (I2.24[) S2 m +i can be an arbitrary 
real number; 

(ii) if £ = 2m + 1 then r(A) = o(l) as X^foo; 

(iii) if £ > 2m + 1 then r(A) has the asymptotic expansion 

(2-28) ^) = -|--|^ + (a^)' A ^°°' 

where the sequence (s, £ — 2m — 2) is determined by the matrix equation 

T{p m +l, ■ ■ ■ ,P0, — £S0j • ■ • , — £Sf_2m-2) T(s m , . . . , Sg) = S m It- m +l. 

Proof. It is clear that the function ip admits the asymptotic expansion (|2.26p if and 
only if c(A) := — X m+1 (p(X) admits the asymptotic expansion of the form (|2.2ip with 
n — £ — m and, moreover, by standard inversion of expansions, this is equivalent for 
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d(X) := l/c(A) to admit the asymptotic expansion of the form (|2.21|) with n 
Now by substituting the expansions for the terms in the formula 



(2.29) 



and applying Lemma 12.111 it is seen that 

/ qr ^ p(A) +er(A) p m p 2m -i+i , ( 1_\ 

(^-«>0) — = p m +i + —: — H ri h o —, , A— »oo, 

v ' A m + r T \ A m \ A / 

where the coefficients are determined by the matrix equation 

(2.31) T(p m+1 , . . . ,p 2m -e+i) T(s m , ...,si) = 

Sm-L£ — m+1 • 

In particular, p(X) appearing in (|2 . 2T[) is a polynomial of degree m + 1 whose 
coefficients satisfy the matrix equality (|2.25j) when £ > 2m. Hence, p(X) in (|2.27p 
can be taken to be the polynomial defined by (|2.24j) in Lemma [2.121 

(i) If £ = 2m then the formula (|2~3"D)) shows that the function r(A) in ([2^27)1 
satisfies r(A) = o(A) as A^>oo, independent from the selection of the real number 
s 2m +i in (I2~m 

(ii) Hi — 2m + 1 then the formula (|2.30[) shows that r(A) = o(l) as A^Poo. 

(iii) If £ > 2m + 1 then the formula ([23D1) shows that the function r(A) in (pT2T)) 
has the asymptotic expansion (I2.28|) with = — ep_,-_i, j = 0, — 2m — 2, 
with coefficients pj as in (|2.31[) . □ 

Associate with the expansion (I2.26|) the (m + 1) x (m + 1) matrix S m as in (|1.4I) . 
Then 



(2.32) S m = 



where T(s, m, 2m) and J m +i are as in (|2.19j) and (|2.20j) . Let the monic polynomial 
p(A) of degree m + 1 be defined by the formula (|2.24|) . where S2m+i is an arbitrary 
real number in the case I — 2m. Then it follows from (|2.25l) that 

and hence 

(2.33) Sm ^^^(l ... X-)slS m 

A — UJ 



so that s m p e N^_(5 m ) (as s m ^ 0), i.e., the negative index of the generalized 
Nevanlinna function s m p(X) is equal to i>-(S m ), the number of negative eigenvalues 
of the matrix S m . 

When the function ip(X) in (|2.26|) is a generalized Nevanlinna function the state- 
ments in the previous lemma can be specified further. The next result shows how 
the classes N K ,_f behave under linear fraction transforms; for this it suffices to 
consider the transform f(X) — > — l/</?(A); as in Lemma 12.131 the result is expressed 
via the function r(A) which will appear in later sections, too. 

Lemma 2.14. Let the notations and assumptions be as in Lemma \2. 13i Then the 
following assertions hold: 
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(i) ip £N K if and only if r G N K _ iy _ (5m ); 

(ii) ip G N Kj _2m if an d only if r G N«-i/_(s m ) <wid T (-M — o(A) as A— »oo; 

(iii) G N Ki _2tn-i if md only if t G N k _„ _(s m ),i T (-M = o(l) as A^?oo; 

(iv) ip G N K _f iwift £ > 2m + I if and only if r G N fc _ ( ,_( l s TO ) ! _™_2 m _2) and 
t(A) = o(l) as A^Foo. 

Proof, (i) The condition 93 G N K is equivalent to — l/<p G N K (^3 5^ 0). Since £ > 2m 
it follows from Lemma \2. 131 that r(A) = o(A), A^oo. Hence the definition of t(A) 
in (|2.27p and [131 Theorem 4.1, Corollary 4.2] imply that —l/(p G N K if and only if 
r(A) is a generalized Nevanlinna function such that its negative index k(t) satisfies 

k = n(—l/(p) — k(sP) + k(t) = v-{S m ) + k{t); 

see (|2.33p . Hence, ip G N K if and only if r G N K _„_(s m ). 

(ii) The statement for I = 2m is obtained now directly from part (i) of Lemma [2.13l 

(iii) & (iv) Let t = 2k or I = 2k — 1 with k > m and rewrite the expansion of 
<p(\) in ((2T261) as follows: 

(2-34) rtA) = -^-...-£^ + £W C(A) = (1), A^oo. 

Then by (|2.13[) 93 G N Kj _ 2 fc+i if and only if C G N K < 1 and, similarly, tp G N K) _2fc 
if and only if C G N re ',o for some ft' < k. Now the expansion in (|2.30j) can be 
rewritten as follows 

rosea p(A) +£T(A) ,Prn, P2(m-k+l) D(X) 

( 2 ' 35 ) 3^1+1 =Pm+l + — + ■■■+ A2fc _ m _! + A2fc _ m _ 1 , A^OO, 

where -D(A) satisfies 



D{X) = C(A) + O f ~ V A^oo; 

s m \XJ 

compare |21[ Lemma 4.1]. The formula ()2.35|) is equivalent to the following expan- 
sion for t(A): 

^) = K¥ + ''' + §E^) + ^tt- 

where eD(X) = EPm+1 C(A) + O (ij as A^oo (and for k = m + 1 the first term 
in the righthand side of (|2 . 36[) is missing). Since here eD(X) G N re »j- for some 
«" G N is equivalent to C G N K / j3 for j = 0,1 (p m +i = 1), the assertion t G 
N K _^_(5 m )._(^_2m-2) for the values I — 2k and I = 2k — \ with fc > m follows from 
(|2T36l) by the inclusions (j243[) . □ 



In the special case ft = and m = the result in Proposition [244] implies [21] 
Theorem 4.2]. Results analogous to that in Proposition 12.141 in the special case 
where I — 2k is even can be found from [T3J Proposition 5.4], [TJl Theorem 5.4]; 
the even case (I = 2k) is easier, since then the expansion (|2.26[) for ip G N K is 
equivalent to <p G N K! _ 2m , see e.g. [14j Corollaries 3.4, 3.5]. 

3. Basic moment and interpolation problems 

In this section solutions of the so-called basic moment and interpolation problems 
will be described. The treatment is divided into two cases: nondegenerate and 
degenerate problems according to det S n 7^ and det S n — 0, where the matrix 
S n is as defined in (|1.4[) . The even and the odd cases of the problems MP K (s,£) 
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and IP K (s,£) will be treated in a parallel way. In the even case {£ = 2n) all the 
statements will be formulated only for the problem MP K (s, 2n), since the problems 
MP K (s,2n) and IP K (s,2n) are equivalent (see Remark I2.9[) . In what follows a 
sequence (s,£) := {si}j =0 is said to be normalized if the first nonzero element of s 
has modulus 1. 

3.1. Nondegenerate basic moment and interpolation problems. Let n = 

[1/2], so that either £ = 2n or £ = 2n + 1. Nondegenerate problems MP K (s,£) and 
IP K (s, £) are said to be basic if the sequence (s, £) is normalized and det Sj = for 
all j < n — 1, or, equivalently, if 

(3.1) s = si = • •• = s„_i = 0, |s„|=l. 

Thus a function ^ € N K is a solution to the nondegenerate basic interpolation 
problem IP K (s,2n), if it satisfies the condition 

(3-2) v(X) = -^-^-----j^ + o^y A^oo. 

If IP K (s,£) is nondegenerate and basic, then the Hankel matrix S n has the form 

/ 

\ s n ... s 2r: 

where all the nonspecified entries are equal to 0. Define the monic polynomial p(A) 
of degree n + 1 by the formula (|2.24j) . where m — n and S2n+i is an arbitrary real 
number in the case of even £. Then it follows from (|2.33p that s n p £ N„_(5 n ). In 
the case of even I = 2n the following result is a corollary of general descriptions 
of A4 K (s,£) given in [18] and [11], and a short proof in this even case £ = 2n has 
been presented in [§]. Here the result both for even and odd £ is an immediate 
consequence of the general transformation result given in Proposition 12.141 

Lemma 3.1. Let (s,£) be a sequence of real numbers such that (|3.1I) holds with 
n = [£/2], let := v-(S n ), let p G N„_ be the polynomial of degree n + 1 defined 
in (|2.24|) , and e = s n . Then MP K (s,£) and IP K [s,£) are solvable if and only if 

(3.3) n>v-. 

If k > V- then the formula 

p(X) +er(A) 

describes the sets M K (s, £) and X K (s,£) as follows: in the even case 

ip e M K (s,£) »t£ N k _ 1/ _ and satisfies (E); 

and in the odd case 

f € M K (s,£) «t£ N K _, y _ i i and satisfies (O); 

if G I K (s,£) »t£ N k _^_ and satisfies (O). 

Proof. Let </? G l K (s,£). Then G N K with the expansion (|3.2p . Since |s n | = 1 
Proposition 12.141 shows that now equivalently 

(3.5) -l/<p(X) = ep(X)+T(\), 
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where r G N K _„_ . In particular, the solvability criterion (|3.3[) and the assertions 
in the even case £ — 2n and the odd case £ = 2n + 1 are obtained from and parts 
(i) and (ii) of Lemma 12.131 respectively, and Lemma 12.141 (i) . 

The last statement for the moment problem MP K (s, £) in the odd case is implied 
by Lemma 12.141 (iii) . □ 

3.2. Degenerate basic problems. Let £ G N and n = [£/2). Degenerate moment 
and interpolation problems MP K (s, £) and IP K (s, £) are said to be basic if det Sj = 
for all j < n and the sequence (s, £) is normalized. Consequently the set of 
degenerate basic moment problems can be divided into two cases as follows: 

(A) Sj = for all j = 0, . . . , 2n. 

(B) There is at least one nonzero moment Sj for some j = n + 1, . . . , 2n. Let m 
be the minimal number for which \s m \ = 1 (n < m < 2n). 

3.2.1. Degenerate basic problems: Case (A). In this case sq = • • • = Sa« = and 
va(S n ) = n + 1. Let us denote vq :— vo(S n ). In the next theorem descriptions of 
the sets of solutions to the degenerate basic problems MP K (s,£) and IP K (s,£) are 
given. 

Lemma 3.2. Let (s,£) be a sequence of real numbers satisfying the assumption (A) 
with n — [£/2]. Then in the even case the problems MP K (s,£) and IP K (s,£) are 
solvable if and only if 

(3.6) either k — 0, or k > vq :— vo(S n ). 

In the odd case MP K (s,£) and IP K (s,£) are solvable if and only if 

(3.7) either k = and sg = 0, or k > vq. 

If k = (and S£ = in the odd case), then it has the unique solution (p(\) = 0. 

If k > vq and v is given by (|2.17[) , then the formula 

(3-8) *») = fg, 

describes the sets M K (s,£) andI K (s,£) as follows: in the even case 

<p € M K (s,£) (p 6 N K _„ and satisfies (E); 

and in the odd case 

tp G ^^^(s,^) sg € N K _ y i i and satisfies (O) ; 

ip G I K (s,f) + € N K _^_ and satisfies (O) . 

Proof. The case k = is trivial, since if <p G No, then so = implies </? = 0. 

Let k > and let ip G N K be a solution of the interpolation problem IP K (s,£). 
Then it follows from (13.21) and (A) that in the even case (£ — 2n) 

(3.9) ^ (A) = (__L_^ a^oo. 

and in the odd case 

(3-10) ^) = - S -^ 2 +o(j± r2 

In both cases 

lim A 2 "+V(A) = 
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and, hence, the multiplicity iToo(<p) is at least n + 1. In view of Remark 12. II one has 

« > TToo(^) > Tl + 1 = V . 

This proves the necessity of the condition k~>vq. Due to Lemma ^.lOl cp admits the 
representation (j3.8[) with G N K _„. In addition, it follows from (|3.9[) and (|3.10[) 
that tp satisfy the assumption (E) if £ is even, and tp + sg satisfy (O) if £ is odd. 

And finally, let £ — 2n + 1 and k > vq = n + 1. Then by Lemma fe.lOl y g N Ki _f 
if and only if £ N K _^i. This completes the proof. □ 

3.2.2. Degenerate basic problems: Case (B). In this case a function tp from N Kr g 
(N K ._£) is a solution to the degenerate basic problem MP K (s,£) (IP K (s,£)), if 

(3-ii) ^) = -^tt & + {jh)> A ^°°' 

with \s m \ = 1 and n < m < 2n, where n — [£/2]. 

Lemma 3.3. Let (s,£) be a sequence of real numbers satisfying the assumption (B) 
with n = [£/2] . Then the problems MP K (s, £) and LP K (s, £) are solvable if and only 
if 

(3.12) k > k :=i/ +v-, u Q := u Q (S n ), V- := v-(S n ). 
Let the sequence (s,£ — 2uq) = {si}fZg 1/ ° be given by the equalities 

(3.13) Sj=Sj +2vo , (j = ) ...,£-2v ) ) 

and let v be defined by (|2.17p . Then the formula (|3 - 8[) establishes a one-to-one 
correspondence between solutions tp of the problem IP K (s,£) and solutions tp of the 
nondegenerate basic problem LP k ,- v {s 1 £—2vq). A similar statement concerning the 
problems MP K (s,£) and MP K - U (s,£ — 2i/q)) is also true. 

Proof. Let tp 6 l K (s,£), so that tp belongs to N K and satisfies p. lip . Now the 
matrix S n takes the form 

0(m — n)x(m — n) 0(m— n) X (2n— m+1) \ 

0(2n — m+1) X (m — n) ^[m — n,n] J 



( 



where 

S[ m -n,n] = ( s i+j)i,j=m-n = 

\ Sm • • ■ ^2r; 

is invertible since s m ^ (\s m \ = 1). It is clear that 

vo(S n ) = m — n > 0. 

To determine the index v_(S n ) consider the following three subcases: 

(Bl) m is even and s m > (denote m = 2k); 
(B2) m is odd (denote m = 2k — 1); 
(B3) m is even and s m < (denote m = 2k — 2). 
Then one can easily check that 

n — k, in case (Bl); 
U-(S n ) = ^ n— k+ 1, in case (B2); 

n — + 2, in case (B3), 
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so that in each of the cases (B1)--(B3) one has 

(3.14) v (S n )+v-(S n ) = k>0. 
It follows from p. lip that in the case (Bl) 

lim A 2fe+ V(A) = -s 2k < 0, lim A 2fe ~V(A) = 0, 

In the cases (B2) and (B3), respectively, one obtains 

lim X 2k+1 ip(X) = oo, lim X 2k ~ l ip(X) = 0; 

A^Voo A^-oo 

lim A 2fe+ V(A) = oo, lim A 2fc_ VM = ~s 2k -2 > 0. 

A^-oo A^>oo 

Hence, in each of these cases, oo is a GZNT of tp(X) of multiplicity ir o{'P) = k; see 
(jUJ. By Theorem [2?2] (or Remark [27T]) this implies the inequality (j3T3fl) . 

Due to Lemma [2.1 01 <p admits the representation (|3.8[) where (p G N K _„. Clearly, 
the expansion (|3.1ip can be rewritten as 

V( X ) = - Xm ^ + l X/-2, + l +° [ X i-2uo + l ) » A ^°°' 

Therefore, </? G X K _ v {s,i — 2v$). These arguments can be reversed to obtain the 
converse statement. 

In the case where £ is odd it follows from Lemma f2. 101 that tp g N Ki _£ if and only 
if tp g N K _ !/) _(^_2 !/0 ). This proves the statement concerning the set A4 K (s,£). □ 

The following Lemma summarizes the results of Lemmas 13.11 and 13.31 

Lemma 3.4. Let (s,£) be a sequence of real numbers satisfying the assumption (B) 
with n = [£/2], let (|3~T2|) hold, let v be defined by (|2J7) . £ — s m , and let 

( Sjyi Sm-\-l \ 

(3.15) ft\)= I det 

aeXa [m-n.n\ S m S m+ i ... S 2 n+1 

\ 1 A ... A 2n+1 J 

where s 2 n+i is an arbitrary real number if £ is even. Then the formula 

e 



(3 - 16) ^ X) = - x^(p(x) + er(x)y 

describes the sets M K (s,£) and X K (s,£) as follows: in the even case 

tp G A4 K (s,£) 4=> tp G ~N K - u _-v and satisfies (E); 

and in the odd case 

G Ad K (s,£) tp + si G ~N K - u _-v,i and satisfies (O) ; 

tp G I K {s,£) <^> tp + S£ G N t _„ _„ and satisfies (O) . 

Proof. Assume that (|3 . 1 2|) holds. In Lemma [3~3l the problem IP K {s,£) was reduced 
to the problem IP K (s,£ — 2vq). By Lemma [3.11 the set I K - V (s,£ — 2v ) can be 
described by the formula 

(3.17) tp.(X) : 



p(X)+er(X) 

where r is a function from the class N re _( !/+!/ _) such that the appropriate condition 
(E) or (O) is satisfied. Substitution of (|57T7|> into flH| yields {335} . 
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Due to Lemma [3?3] and Lemma IQI o? belongs to A4 K (s, t) with odd £, if and only 
if the function r belongs to N K _ („+!,_) and satisfies the condition (O). □ 

4. SCHUR ALGORITHM 

The present approach to the degenerate moment problem is based on the follow- 
ing reduction algorithm which for the nondegenerate case with even index £ was 
considered in [S]. 

4.1. One step reduction for moment problems which are not basic. Let 

s = { s jYj=o b° an arbitrary normalized sequence of real numbers and let S n — 
(si + j)f -_ be the Hankel matrix as defined in (|1.4p . Assume that S n ^ and 
consider a sequence of normal indices of S n , 

(4.1) 0<n 1 <---<n N <n + l 
which are characterized by the conditions 

(4.2) detS„ 3 _i^0 (.7 = 1,..., iV). 

In particular, the first normal index m is the minimal natural number such that 
det5 ni -i 7^ 0, or , equivalently, that 

(4.3) s = si = • • • = s ni _2 = 0, s ni _i^0. 

Note that the first normal index satisfies ni(= un) = n + 1 precisely when the 
moment problem is nondegenerate and basic and that there are no normal indices 
for moment problems which are degenerate and basic; see Section [3] 

In this section it is assumed that the moment problem is not basic, i.e., one has 
n\ < n. Let the sequence (s, £) = {Sj}j_ be normalized and denote 

£i = sgns ni _i = ±1. 

In this case a function tp £ N K ._£ is a solution to the moment problem A4 K (s,£) if 

( 4 - 4 ) ^ A )=-^ir-^T+T X^ +0 {x^)> A ^°°' 

Then —1/ip G N K and, moreover, by part (in) of Lemma 12.131 — 1/V admits the 
representation 

(4.5) - 1/V(A) = eipi(X) + a^i(A), 

where Px(X) = Pn} A™ 1 + • • ■ + p^ is a monic polynomial of degree n\ (pni = l)j 
defined by the equation (12.241) with m = n\ — 1, and ai(> 0) is chosen in such a way 
that the sequence (sW, £ — 2ni) = (s^^Zq™ 1 defined by the expansion of ^i(A) 

is normalized. Moreover, by Proposition 12.141 (iii) ip\ is a generalized Nevanlinna 
function from the class 

N re _ Kl) _(^_2n 1 ) ) K i ■— V-{S ni -l)- 

As was shown in Lemma [2.13l the moment sequence (s^, £~ 2n\) is uniquely defined 
by the matrix equations 

(4.7) T(s ni _i, . . . , Sj +2ni )T(p ( V, . . . ,p§\ -eia\s^\ . . . , -ei_a\sf>) = ei/j +ni+2 , 
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where < j < £ — 2n\ . 

The above considerations yield the following result. 

Proposition 4.1. Let S n be a Hankel matrix, let n\ be the first normal index of 
Sn, n-i < n, let the monic polynomial pi(X) — Pni X ni + • • • + Pq and the induced 
moment sequence (s^,£-2n 1 ) be defined by e x — s rii — 1j K l ■ — ^—{Sni—l)- 

Then the formula 

(4-8) ip(X) = T x [ Vl (A)] := ~ £l 2 

establishes a one-to-one correspondence between the sets A4 K (s, £) and M K _ Kl (s*- 1 ), £- 
2ni) as well as between the sets I K (s,£) and I K - Kl (s' 1 ), £ — 2rii). 

The normal indices of the induced Hankel matrix S„l ni = (s^-)"^o can be 
derived from the normal indices of the original Hankel matrix S n . This is given in 
the next Proposition. 

Proposition 4.2. Let n\ < n2 < ■ ■ ■ < n^(< N + 1) be all normal indices of the 
Hankel matrix S n . Then the normal indices of the induced Hankel matrix <5„_ ni 
are 

n% — tii < ■ ■ ■ < n N — n\. 

Proof. It follows from (|A.3|) in Lemma IA.3I that det $£_ ^ if and only if i — 
U2, . . . , rt 7y (n\ < i < n). □ 

Now applying Proposition 14.11 to the matrix S^l n one construct a polynomial 

P2 and a Hankel matrix S n _ n2 . After N inductive steps one obtains the Hankel 

matrix S^} n of induced moments and subsequent application of Lemma lA.3l yields 
the following 

Corollary 4.3. Let S^_ n . be the Hankel matrix of induced moments after j steps 
(1 < j < A^). Then the set of normal indices of the Hankel matrix S^}_ n . takes the 
form {nk — n-j}k=j+f Moreover, for all i such that nj < i < n one has 

(4-9) v±(sV\ j )=v ± (S i )-v ± (S nj - 1 y, 
(4-10) MS^ ni ) = MSi)- 

In particular, the matrix S£j n has no normal indices anymore, that is det — 
for all i such that hn < i < n. 

Proof. The first statement is implied by the formula ()A.3|) in Lemma IA.31 The 
formula (|4.9p can be obtained by induction. Indeed, for j = 1 the statement is 
contained in (|A.2[) of Lemma [A.31 Assume that (|4.9p holds for some j (1 < j < N) 
and all i such that nj < i < n. Then it follows from (|A.2|) that 

In view of the induction assumption this yields 

= v±(Si) - v±{S nj+1 -x). 
The formula (|4. 10[) is immediate from (|A.3I) in Lemma [A. 31 □ 
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4.2. Algorithm. Let us define a sequence Ki < • • • < kn by the equalities 

(4.11) K S = v-{S nj -x), j = l,...,N. 

Due to Proposition 14.11 on each step one obtains a linear fractional transformation 

(4.12) w-'W = 7i[v,Wl== ftW+ ^, wW . 

where 

e j = sign si7_ 1 2._ 1 _ 1 (= ±1), aj > {0 < j < N — 1). 

The transformation 7} establishes a one-to-one correspondence between the sets 
€ A^k-^.^sW- 1 ),^- 2nj_i) and Vi e M K _« 3 .(s«,^ - 2n,-). Let W 3 {X) be 
the matrix 

(4-13) W j (X)= ( ° , j EN. 

\ £ J a 3 a 3 ) 

associated with the transformation Tj (1 < j < N). 

After the j-th step we obtain the following representation for the solution ip of 
the moment problem MP K (s,£) 

ip(X) =Tio---oTj[<pj(X)\ 



(4.14) ei ei£ 2 a? ^'-l^ 



Pi (A) - p 2 (A) - ••• - p j (X) + e j aj<p j (X) 

where the last formula stands for the continuous fraction expansion (this shorthand 
notation is often used in the literature). The resulting matrix W\\ ,ji(A) of the linear 
fractional transformation in (14.14[) coincides with the product of the matrices Wi 
(1 < * < j) 

(4.15) W [1>j] (X) = W 1 (X)...W j (X) (j<N). 

Theorem 4.4. Let n\ < ••• < n^{< n) be a sequence of all normal indices of S n 

fw^(X) w^(X)\ 

and let the matrix ,-ifA) = LV )% \ be qiven by (14.151) . Then for 

[MV \ w 2i (A) wg\X)J 
every j < N — 1 i/ie formula 

,„m r umi- 4? (A V, (A) + w ^ (A) 

(4.16) V(A) = % ]( A)^(A)j .= { ) — — 7-y— , 

<i {X)w(X)+wg>(X) 

establishes a one-to-one correspondence between the sets X K (s,£) andI K - Kj (s^\ £ — 
2rij), where is defined recursively by (14.71) and Kj = v-{S nj -\). Moreover, 

ip e M K (s,£) ^ ipj G A4 K _ Kj (s w ,.g- 2nj) 

In the case, when det S n — the statement remains valid for j = N . 

Proof. The proof is obtained by successive application of the Schur algorithm de- 
scribed above and Propositions l4. 11 14.21 and Corollary |4.3l to the problem MP K (s, £). 
In the nondegenerate case this process terminates when j = N — 1, since n < n^. 
In the degenerate case Propositions 14.11 can be applied one more time, since n > 
npf. □ 
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To find an explicit form of the matrix W^ri.flCA) let us define the so-called poly- 
nomials Pj(X) and Qj{\) of the first and the second kind, respectively, as solutions 
of the difference equation 



(4.17) £j-iejaj-iUj-2 - Pj{\)u.j-\ + a-jUj = ( j = 1, N), 

with the initial conditions 

P (A) = 1 ! P!(A)- Pl(A) 



(4.18) 

Qo(A) = 0, <3i(A) 
As is easily seen from (|4.17l) 



fli 
£i 
ai 



j 3-1 

de s p j = Ui > de s Qj = ni ^ - 1S> - 

»=1 j=l 

Theorem 4.5. The resolvent matrix Wuji(\) in (|4.15l) admits the following rep- 
resentation 

- £ 3 a jQj-iW -Qj(A) N 

where Pj and Qj (1 < j < N) are polynomials of the first and the second kind 
associated with the matrix S n via f|4. 1 T[) . (I4.18p . 

Proof. For j = 1 the formula (|4.19l) coincides with (14. 13[) . Proceed by induction 
and assume that (|4.19j) holds for j := j - 1. Then it follows from (|4.19[) . (|4.13[) and 
the difference equation (|4. that 



(4.19) WWA) 



W[i,,](A)( J )=%-i](WW( ! 



Ej_iOj-iQj-_2(A) -Qj_i(A)\ / -Ej/a,j 
- i _ia i _ 1 P i _ 2 (A) -Pj-i(A) / ^ p 1 {X)/a ] 

1_ ( ej-i£ja,j-iQj-2(\) - Pj{\)Q 3 -i(XT 
a 7 I, -£ j -ie j a j - 1 P j - 2 (X) +PjPj-i(X) 



3 

Due to the difference equation (|4.17l) 

o ^ f-Qj(X) 



(4-20) 



Hence one obtains 



(4.21) W[i,i](A)( J ) =W M _ 1] (A) 



1 \._w. . ..mf o \ /^-Qi-i(A) N 

The formulas (|i^D|) - (|4T5T]l prove (|4~lH . □ 



Remark 4.6. The recursion algorithm for Nevanlinna functions is well known (see for 
example [T] ) . The formula (|4.19l) for the resolvent matrix can be found in [53] , were 
truncated moment problems were studied. The operator approach to such problems 
was presented in |15) . |16) . In the indefinite case this algorithm was studied by M. 
Derevyagin in 9 a , formulas (|4.19j) for the matrix W(X) and the statement of the 
theorem for the nondegenerate even moment problem were proven in I .if) . The 
linear fractional transformations similar to Tj (so-called Schur transform) has been 
studied by D. Alpay, A. Dijksma and H. Langer in [3J, [4]. 
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5. Description of solutions 

In this section we find a solvability criterion and describe the set of solutions 
of the problems MP K (s,£) and IP K (s,£) in the general setting. As well as in the 
case of basic problem we will distinguish non-degenerate problems and two types 
of degenerate problems: 

(A) rankS„ = n N = rankS njv _i; 

(B) rankSVj > njy = rankS nN —i. 

5.1. Non-degenerate moment problem. 

Theorem 5.1. Let (s,£) be a sequence of real numbers such that detS n ^ 0, Then 
the moment problem MP K (s,£) is solvable if and only if 

K>V-. 

The sets A4 K (s,£) andT K (s,£) are parametrized by the formula 

i»2i (A)r(A) + w 22 '{X) 

where in the even case 

if G M K (s,£) »t£ N k _„_ and satisfies (E) ; 
and in the odd case 

<f G I K (s,£) <St£ N k _ 1/ _ and satisfies (O) ; 
if G M K (s,£) <4> t G N K _ l/ _4 and satisfies (0) . 
Proof. The proof is obtained by compilation of Theorem 14.41 Theorem 14.51 and 



Proposition 13. II □ 

This result seems to be new even for the odd Hamburger moment problem. 

Corollary 5.2. Let s , si, • ■ • , S2n+i oe r ^ numbers, such that S n > and det S n ^ 
0. Then the moment problem (|l.ip with £ = 2n + l is solvable and the formula (ll.5[) 
describes the set of solutions of when r is ranging over the class No,! and 

satisfies the condition (O). Moreover, tp G To(s,2n + 1) if and only if r belongs to 
No and satisfies (0). 

The following example shows the importance of the condition r G No.i in Corol- 
lary EH 



Example. Let so = l,sj = 0. Then p(A) = A and the set of solutions of the problem 
Xo(s, 1) is described by 

p(A) +t(A) 

where r G No and r(A) = o(l) as AAoo. The function r(A) = — - — t=- belongs 

to the class No and satisfies the condition r(A) = o(l) as A^>oo. Therefore, 

-On(l + \/A) 

¥HA) = f= 

iAln(l + y/X) - 1 

is a solution of the problem /-Po( s , 1) 

-1 / 1 

A- + ° A 2 



</?(A) = — + o ( — ) as AAoo. 
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However, t ^ No,i (see [20]) and, hence, (p is not a solution of the moment prob- 
lem ([TTT]) . 

5.2. Degenerate moment problem. Case (A). In Theorem 15.81 solvability cri- 
teria for degenerate moment problems with minimal negative signature n = V—(S n ) 
are given. We start with two lemmas. 

Lemma 5.3. Let (s, 2n) be a sequence of real numbers such that det S n = 0. Then 
the Hankel rank Un of the sequence (s, 2n) coincides with the largest normal index 
un of the Hankel matrix S n . 

Proof. By Frobenius Theorem (see [IH Lemma X.10.1]), if r is the smallest integer 
t (0 < r < n), such that (|1.7I) holds, then det <SV-i 7^ 0. Hence r is the normal 
index of S n . Moreover, r is the largest normal index of S ni since the vectors 
(sj, . . . ,Sj +n ) T , (0 < j < n N ) in (|1.7p are linearly dependent. This implies that 
det S n = for all j > r. □ 

Lemma 5.4. Let (s, 2n) be a sequence of real numbers such that detSVi = 0, let 
71 TV be the largest normal index of the Hankel matrix S n and let S n admit a Han- 
kel extension S n+ i, such that v_(S n +i) = i>-(S n ). Then there are real numbers 
cyq, . . . , a nN -i, such that 

(5.2) Sj = a>oSj^ nN H ha njv _is 3 _i (n N < j < 2n + 1); 

(5.3) s 2n > Oi s 2n 

Proof. Let us set Vj — (sj, . . . , s J+ „) T , (0 < j < n + 1). Since V-(S n ) = v^(S n +\), 
it follows from Lemma IA.2I that 

v n +i € span (vo, . . . , v n ), 

i.e. there is c e C n+1 such that 

(5.4) v n+1 = S n c. 

By Lemma 15.31 there are real numbers ap, . . . , a r _i, such that 



(5.5) 



•so 



where a := 



S2n/ 



a r -i 

\0( n -r+l)xl/ 



This together with (|5.4p implies, in particular, that 



(5.6) 



Sr-\-n-\-l 



C V r 



c* S n a 



S2n+1 ) a. 
\n+l 



Denote by V the (n+ 1) X (n + 1) forward shift matrix V = (^,j+i)™j=i ■ Then (|5?2 
and (|5.6p imply 



/ s r+l \ 


■( 


\S r+n+ iJ 





•Si 



»n+l 



S2n+1 , 



Iterating these calculations one obtains for (0 <)j<n 



(5.7) 



Vj+l 



( 31 

\Sn+l 



S2n+1 , 



-r+l 



a, 
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which proves (|5.2j) . Setting in (|5.7|) j = n one obtains 

v n+1 = S n V n - r+1 a. 
Then it follows from Lemma [A. 21 that 

S2n+2 > (V a) b n V a 

(Sn+1 
: 
S2n+1 , 

= «0 s 2ri-r+2 H + ar-lS2ri+l- 

Hence (|5.3p holds and this completes the proof. □ 
This motivates the following definition which in the definite case was used in [8] . 

Definition 5.5. A sequence (s,£) — {sj}j =Q with the Hankel rank r = rank s is 
called recursively generated, if there exist numbers ao, . . . ,a r _i, such that 

(5.8) Sj = a Q s 3 - r H h ctr-iSj-i (r < j < £). 

Theorem 5.6. Let (s,£) be a sequence of real numbers such that detiSVi = 0, 
n — [1/2], let ni < ■ ■ ■ < un be all normal indices of the degenerate Hankel matrix 
S n , let (s( N ',£ — 2njy) be a sequence of induced moments determined by successive 
application of (|4.7[) . and let k = v_(S n ), = v-(S„ N -i). In the case when 
I = 2n is even the following statements are equivalent: 



(i 
(ii 
(hi 
(iv 

(v 
(vi 



The moment problem J\4 K (s,£) is solvable; 

The moment problem Aia(s^ N \£ — 2n^) is solvable; 



s {N) - ■ ■ ■ - s (N) -0- 



S n admits a Hankel extension S n +i such that ^(SVi+i) = V—(S n ); 
(s, £) is recursively generated; 
rank S n — n^; 
If £ — 2n + 1 is odd, then 

(i) (ii) 4^ (Hi) 4^ (iv') ^ (v) (vi) 

where ( iv ') and ( vi ') take the form: 

(iv') there exists a real number S2n+2, such that v>-(S n +i) — v-(S n ); 

(vi') rank S n = njy and s^_l n = 
If one of the above conditions holds, then n = kjv and the moment problem A4 K (s, I) 
has the unique solution given by 

Qn(X) 



(5.9) ip(X) 



Pn(X) 



Proof. Even case, (i) =>• (iii) If the moment problem MP K (s, 2n) is solvable then 
by Theorem 14.41 the problem MP K _ KN (s( N \ 2(n — n^)) is also solvable. If (iii) is 
not in force then it follows from Proposition 13.31 that 



In view of Corollary 
(5.10) M^L) = v-(S„) - V-(S nN -i) = k-k n . 



24 



V.A. DERKACH, S. HASSI, AND H.S.V. DE SNOO 



Therefore, 



-(CL)>-(^L) + ^o(5i w L), 



which implies vo(S^} n ) = 0. But by the same Corollary 14.31 vo(S^J ns ) = 
vo{S n ) 0. 

(iii) =>> (ii) If (iii) holds then V-(s£) nif ) = and by Theorem O^jv (A) = is 
the unique solution of the problem M P (s( N \ 2(n — n^)). The equality K = is 
implied by (|5~TU1) . 

(ii) => (i) This follows from Theorem |4~4"1 

(iii) => (iv) If (iii) holds then in view of Theorem l4.4l and Theorem 14 . 5 1 the unique 
solution of the problem MP K (s, 2n) is given by 

Since <p is rational of rank ip = it admits the asymptotic expansion (|1 .3[) for 
every n, in particular, there exist S2n+i, S2n+2 € K such that 

1 



^ (A) = A A^ 



- s 2n+2 
\2n+3 



\2n+3 



(A— >oo) 



and the corresponding Hankel matrix S n +i has the same rank as S n . This implies 
that i/_(5 n+ i) = y-(S n ). 

(iv) =>• (v) Assume that S n admits a Hankel extension S n +i such that i/_ (5„+i) = 
i/-(S n ). Then by Lemma T5.3I there are cto, . . . , a nN -%, such that 

(5.11) Sj = a Sj-njv H 1" "n^-iSj-i ("JV < j < 2n + 1). 

By Definition 15.51 this means that (s, 2n) is recursively generated. 

(v) =>■ (vi) The condition (v) implies that 



(5.12) 



£ span 



so 



»T1JV— 1 



for all j < n, and therefore, rank = n^. 

(vi) =>■ (iii) If rank S n = tin then Vo(S n ) = n — + 1. By Lemma I A. 31 



and, hence, = for all j such that < j < 2(n — tin). 

Odd case. In the odd case the proof of the equivalences (i) (ii) (iii) =>■ 
{iv') is pretty much the same. 

(iv') => (v). Assume that there exists S2 n +2 such that z/_(5„+i) = V-(S n ). 
Then by Lemma f5T3l there are ao, . . . , a nN -\, such that (|5.1ip holds for all for 
all (njv <)j < 2n + 1. This implies, that the sequence (s, 2n + 1) is recursively 
generated. 

(v) =>■ (vi') The statement (v) implies that (|5.12l) holds for all j < n + 1 and 
hence there exist /3i , . . . , j3 nN , such that 




' Slip, 



TRUNCATED MOMENT PROBLEMS 



25 



Therefore, rank S n = tin. Let us set 

S2n+2 = P\S n +l + ' ' ' + p nN S n + nN . 

Then rank 5 rl +i = rank S n = and by (iii) we obtain 

„W _ _ JN) _ (TV) _ „ 

°0 ,, 2n+l-2»ijv — a 2ri+2-2njv — 

(vi') => (iii) If rank S n = njq then it was shown above that Sj = for all 
j < 2(n — npf). Now (iii) holds since s 2 ^ 1 _ 2nN = 0. □ 

Remark 5.7. If p is a rational function of degree r and p has the asymptotic 
expansion (|1.3[) . then for n > r— 1, by Kronecker theorem rank S n = r and <p € N K , 
where k = v_(S n ) (see [T^l Theorem 16.11.9]). Then by Theorem 15.61 the problem 
IP K (s, 2n) has a unique solution. 
Now, let ip € N K be such that 

tp(X) = <p(X) + o f^2f+T ) as A ^°°- 

Then both 93 and -0 are solutions of the problem IP K (s, 2n) and hence ip(X) = ip(X). 
This proves the rigidity result for generalized Nevanlinna functions obtained in [5] 
and proved originally by Burns and Krantz for functions from the Schur class, [7]. 

In the next theorem we describe solutions of the problems MP K (s, £) and MI K (s, t) 
in the case where the rank of the Hankel matrix S n coincides with the Hankel rank 
of the sequence. 

Theorem 5.8. Let (s,£) be a sequence of real numbers such that det S n = (n = 
[1/2]) and let 

rank S n = 11 ^ . 

Then the problems MP K (s 1 £) and MI K (s,£) are solvable if and only if: 

(i) either k = V—{S n ) and the equivalent conditions of Theorem \5.6\ are satis- 
fied; 

(ii) or k > v-(S n ) +vo(S n ). 

If k> v-(S n ) + va{S n ), W[i,n^(X) is given by (|4.19jl . and 

(5-13) W(X) = W [liN] (X) ( I x i Q ) 

then the formula 

(5-14) <pW=T wW [t(\)], 

describes the sets M K (s,£) and X K (s, £) as follows: in the even case 

p e M K (s,£) «t£ N k _„ and satisfies (E) ; 
and in the odd case 

p € 2 K (s,£) t + s^L € N K _„ and satisfies (0) ; 

if G M. K [s,i) r + s^\ nN £ N K _ I/j i and satisfies (0) . 
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Proof. It follows from Theorem 14.41 that the problem MP K (s,£) is solvable if and 
only if the basic problem MP K _ L/ _(§ n _ 1 ^(s( N \ £—2n^) is solvable. Since rank S n = 
rankSjj^-i, then u±{S n ) = v±(Sn N -i)i ari d one obtains from Corollary I4.3I that 
rankjS^-j = 0. Hence the basic problem M. R - V _ (s^ N \2(n — njv)) is of type (A). 
The rest of the statements are immediate from Lemma I3.2I □ 

5.3. Degenerate moment problem. Case (B). In this subsection we give solv- 
ability criteria and describe solutions of the problems MP K (s,£) and MI K (s,£) in 
the case where the rank of the Hankel matrix S n is greater then the Hankel rank 
of the sequence (s,£). 

Theorem 5.9. Let (&,£) be a sequence of real numbers such that det S n = (n = 
[£/2]) and let 

rankS'n > riff, 

Then the moment problem A4 K (s,£) is solvable if and only if (jl . 13[) holds. 

Let V- := V-(S n ), vq := fo(S n ), let v be defined by (|2.17j) . let (s^ N \£— 2un) be 
the sequence of induced moments after N steps, let the integer m be defined by 

sf } =0forj<m; e := »W + 0, 

let the polynomial p be defined by the formulas Q3.13P and (|3.15p , where n := n—npf, 
let the matrix-valued function W[i^n](\) be given by (|4.19p . and finally let 

If k satisfies (jl . 13|) . then the formula (|5.14l) . describes the sets A4 K (s, £) andI K (s,£) 
as follows: in the even case 

(f £ M. K (s,£) Ore N K _„_j,_ and satisfies (E) . 

and in the odd case 

ip € I K (s,£) r + SfZlnx ^ N K - U - U _ and satisfies (O) ; 

<f G M K (s,£) t + sf_l G N s _„_„_,i and satisfies (0) . 

Proof. It follows from Theorem I4.4I that the problem MP K (s,£) is solvable if and 
only if the basic problem MP K _ u _(s n _ 1 )(&^ N \£ — 2un) of type (B) is solvable. 
By Proposition I3.3I the latter problem is solvable if and only if 

Due to Corollary I4.3I 

(S i n N } nN ) = MSn), *>- (S { n N l lN ) = U-(S n )-U- (S„ N - 1 ) . 

These proves the criterion (| 1 . 13[) . The statement for the problem IP K (s, £) is proved 
by the same reasonings. 

The second part of the proof is implied by Theorem 14.41 Theorem 14.51 and The- 
orem 13.41 □ 
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Appendix A. Some lemmas on matrices 

A.l. Some lemmas for block matrices. The verification of the first lemma is 
left to the reader. 

Lemma A.l. The inverse of the invertible block matrix A of the form 





( ° 





.4 


H 





A 2 2 


A 






A* 

^23 


A: 




~^13*(^33 — ^23^22 1 ^2 3 )^i 3 1 —A^A^A^ A 13 

-A^A 23 A^ 
Aii 

The following lemma extends a well-known result for nonnegative block matrices. 
Lemma A. 2. Let A be a Hermitian matrix of the form 

A=(g* c) ' AeC nxn , C eC mxm 

such that i/- (A) — f_(A). Then 

(i) ran£? C ranA; 

(ii) if Bh = Ag for some h € C"\ g € C n , then h*Ch > g*Ag. 

Proof. First assume that A is invertible. Then the identity 

A B\_( I 0\ (A \ (I A- X B 

B* C')~\B*A- 1 II [o C-B*A- 1 b)\Q I 



A = 



and the assumption v-{A) = v-(A) shows that C — B* A 1 B > 0. If A is not 
invertible, then there is a block decomposition 



0/ ' \B h 

where Ao is invertible. Note that v~{A) = v-(Ao). Now the previous statement 
can be applied to the block matrix 

'A B S 
B x 

K Bl B\ C t 

with the invertible matrix Aq in the left upper corner, since v~{A) = i/_(Ao). Thus, 
one concludes that 



( 













^(0 = c _ B f A ^ Bo )>0- 



This implies B\ = and 

(A.l) C-B*A^B >0. 

The identity B\ = means that ran B C ran Ao = ran A proving the range inclusion 
in (i). On the other hand, it follows from (|A.1|) that the vectors h and go for which 
Boh = Ao.go satisfy the inequality 

h*Ch > h* BqAq 1 B h = g* a Aogo = g* Ag. 
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Now (ii) is implied by this inequality. 

A. 2. Some results for Hankel matrices. 



□ 



Lemma A. 3. Let S n be a Hankel matrix, let ri\ be the first normal index of S n , and 
let the polynomial p i (X) = pnj \ ni + ■ ■ '+Pq and the moment sequence (s^,^— 2ni) 



be defined by (|4.7p . Then for all i = 0, n — n± 

(A.2) v±{S\ 1] ) = v±(S i+nx ) - v±(S ni -i); 

(A.3) voisV) = v (S i+ni ). 

Proof. Consider the equation (|4.7p with j = 2i. Multiplying it both from the left 
and from the right by the matrix one obtains the equality 



(A.4) 
where 



A t2 



12 



Bn B 
B{ 2 B 2 2 



eil. 



2i+ni+2, 



Al2 = 



Sni-1 \ 



£(i+l)x(ni+i+l) 



(A.5) 



... s ni — \ ... s ni -\-i—\ J 

A 2 2 = S ni + l eC («i+i+l)x(n 1 +i+l) j 
-oil — — S\a 1 Ji + ib i Ji + i G IL V 

and B 12 , B 22 are some matrices from c( l+1 2 x(i+ ™ 1+1) and c( i+ni+1 ) x ( i +" 1+1 ), re- 
spectively. Let us decompose the matrices A\%, A22 as follows 

li 2 = (0 (m)xni A 13 ), A 13 = .■ : eCt'+^t'+D, 

— 1 ■ ■ • Sn\-\-i— 1 / 

'A 22 A 7A \_(S ni -i A 23 \ A 23 G C" lX ( J+1 ), 



A 22 = 



A 33 



4' ! 



-4, 



Ao 3 G C^+Dx^ 1 ). 



*23 ^33, 

Since the matrices A13 = A* 3 , A22 = S ni —i are invertible, it follows from Lemma 
Al that the matrix 





f 





.4 







A 22 


.4 




\A 13 


^23 


-4 



is invertible, and the left upper corner Bn of the inverse matrix A 1 is given by 
(A.6) B n = (A~ 1 ) n = -A^(A 33 ~ A* 23 A^A 23 )A^ 

It follows from (fO|) - (fO)) that 

(A.7) alJ l+1 s\ 1] J l+1 = -£iB n = -(I" 1 )!! = ^(Ass - A^ 1 ^)^ 1 
and hence 

(A.8) v±(.S?>) = v±(A 33 - A* 3 A 22 1 A 23 y, 
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2!) 



(A.9) 



v {s[ l) ) = v {A 33 -A* 23 A 22 1 A 23 ). 



The numbers of positive (negative, zero) eigenvalues of the Schur complement A 33 - 
A 23 A 22 A 23 can be calculated by the formulas 

(A.IO) ^±(A 33 - A* 23 A 2 ^A 23 ) = v±{A 22 ) - v ± {A 22 )- 

(A.ll) ^0(^33 ~ A* 23 A 22 1 A 23 ) = MA22) - v$(A 22 ). 



Since A 22 — Si+ ni and A 22 = S ni -i is invertible the statements of the Lemma lA. 31 
arc implied by (|A^8 ]l - (|ATTT]l . □ 
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